Yang-Lee zeroes for an urn model for the separation of sand 
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We apply the Yang-Lee theory of phase transitions to an urn model of separation of sand. The 
effective partition function of this nonequilibrium system can be expressed as a polynomial of the 
size-dependent effective fugacity z. Numerical calculations show that in the thermodynamic limit, 
the zeros of the effective partition function are located on the unit circle in the complex z-plane. 
In the complex plane of the actual control parameter certain roots converge to the transition point 
of the model. Thus the Yang-Lee theory can be applied to a wider class of nonequilibrium systems 
than those considered previously. 
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/. Introduction. — Nonequilibrium phase transitions 
have recently attracted increasing theoretical interests. 
One of the motivations is the idea that some of the con- 
cepts of equilibrium phase transitions might apply to 
nonequilibrium situations as well. In particular, there 
have been recent investigations Q in the applicability of 
the Yang-Lee theory to nonequilibrium systems. Usually, 
those systems are modeled in terms of a master equation. 
Most of the physics is contained in the prescribed tran- 
sition rates. When varying the control parameters, the 
system may exhibit a nonequilibrium phase transition in 
its stationary state. There are very few models for which 
such transitions can be described analytically. Moreover, 
it is not obvious to find a nonequilibrium function play- 
ing a similar role to the equilibrium partition function 
and for which the Yang-Lee strategy might be extended. 
Two families of dynamic models have been discussed so 
far in the literature 0. 

(i) Most of the works concern driven diffusive sys- 
tems 0, S These models have a unique steady state, 
for which one can compute exactly the probability dis- 
tribution of microstates. The normalization factor, de- 
fined as the sum over these probabilities, is shown to 
play the same role as the partition function in equilib- 
rium systems. It is a polynomial in the control param- 
eter, the latter being a constant transition rate, i.e., a 
size-independent quantity. The roots of this polynomial 
in the complex plane of the control parameter accumu- 
late, in the thermodynamic limit, at the transition point. 
As in the case of equilibrium transitions, the way these 
zeroes accumulate nearby the real axis informs about the 
type of transition (first order if the line of zeroes is per- 
pendicular to the real axis, and second order if it forms 
an angle of 7r/4 with this axis). 

(ii) Other studies have been performed for directed 
percolation models on a square lattice in (1 + 1) dimen- 
sions It was suggested that survival probability P(t) 
for a finite propagation time t might play the role of a 
partition function in this case, with the occupation prob- 
ability p as the complex control parameter. However, the 
survival probability, that is actually the order parameter 
in this problem, does not have the "standard" properties 



of an equilibrium partition function [l(. Therefore it is 
not surprising that the distribution of its complex zeroes 
has a more complicated structure. It was also shown that 
the minimal distance d(t) between the complex zeroes of 
the survival probability and the critical point p c scales 
for large t as t -1 /"!! , where v\\ is the critical exponent for 
the temporal correlation length. 

The above examples suggest that it is rather straight- 
forward to apply Yang-Lee theory to nonequilibrium 
phase transitions, at least as long as one can construct 
the appropriate effective partition function (EPF). How- 
ever, an implicit assumption of the above approaches is 
that the EPF can be expressed as a polynomial of a size- 
independent control parameter. The purpose of our letter 
is to show, on a simple model, that this should not be nec- 
essarily the case. Indeed, our model has an EPF that has 
a polynomial structure in terms of a size- dependent effec- 
tive fugacity, and not in terms of some size-independent 
parameter. The applicability of the Yang-Lee approach 
to such a case is thus highly questionable. We show that 
the Yang-Lee strategy for this model still works and thus 
can be generalized to a wider class of nonequilibrium 
phase transitions. 

This Letter is organized as follows. In Sect. II the 
model is defined and the EPF is introduced as a polyno- 
mial of an effective size-dependent fugacity. In Sect. Ill, 
the zeroes of the EPF are numerically studied. They 
are shown to form a much more complicated structure 
in the plane of the complex fugacity than in the above 
discussed examples (i) and (ii). Nevertheless, it is still 
possible to extract some information about the nature of 
the phase transition in our model by investigating the 
zeroes in the complex plane of the control parameter. 
In particular, a subset of zeroes converges to the tran- 
sition point. Moreover, the convergence seems to obey 
standard finite-size dependence arguments. In the Con- 
clusions (Sect. IV), we argue that the behavior of zeroes 
in our model is an indication of difficulties that might be 
typical to more general nonequilibrium systems, in par- 
ticular those where transition rates between microscopic 
configurations are not constants (i.e., size-independent 
quantities), but state dependent. 
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77. The model. — The model we consider was intro- 
duced to describe spatial separation of vibrated sand 
and is a generalisation of Ehrenfest's urn model 0. In 
this model TV particles are distributed between two urns, 
the first urn containing M particles and the second one 
N — M. The dynamics is defined as follows. At each 
time step, one particle is chosen at random in one of 
the urns. Then, with a probability that depends on the 
number of particles present in the chosen urn, i.e., with 
a state- dependent transition rate, this particle moves to 
the other urn. Correspondingly, the flux F(n) of parti- 
cles leaving a given urn at a certain time depends on the 
fraction n of the total number of particles in the given 
urn at that moment. This model is thus by construction 
mean-field like. The master equation for the probability 
distribution p(M, t) that there are M particles in a given 
urn at time t writes (J: 
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Its stationary solution is found to be 



N 



Z "l\ Fit) 



p(M,t). (1) 



(2) 



where the normalization factor is 
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In the spirit of models of type (i), we shall assume that 
Zm plays the role of an EPF. 

This model describes the transition between a sym- 
metric distribution of the particles in the two urns, as- 
sociated with a single peak of the probability distribu- 
tion at M = TV/2 (for TV even), to a symmetry breaking 
state described by a bimodal distribution with peaks at 
M = N(l/2 ± e). The order parameter e measures the 
difference in the occupancy of the two urns. To produce 
this symmetry breaking it suffices that the flux function 
F(n) has a single hump [y, |9j- The simplest possible 
choice for F(n) having this property is 



F(n) — n exp (—An) , 



(4) 



which corresponds to a state- dependent transition rate 
exp(— An). Thus the problem is characterized by a single 
control parameter A. In the thermodynamic limit TV — > 
oo, this symmetry breaking corresponds to a second- 
order nonequilibrium phase transition. In this limit, 
the probability distribution becomes 5-peaked around the 
macroscopic stable state, that is determined by the condi- 
tion that the flux of particles directed from the first urn to 
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FIG. 1: Zeroes of the EPF JSJ in the complex z-plane for 
three values of TV. The continuous line is the unit circle. The 
inset illustrates the behavior in the vicinity of z = 1. 



the second one equals the flux of particles from the second 
urn towards the first one, F (1/2 - e) = F (1/2 + e). It 
follows that in the thermodynamic limit for A < A c = 2 
the stationary state is the symmetric one, while for 
A > A c = 2 the equipartition of particles is broken, i.e., a 
second order nonequilibrium phase transition takes place 
at ^4 = A.Q = 2 . 

777. Analysis of the zeroes of the EPF. — With such a 
choice of the flux F(n), one may rewrite the normaliza- 
tion factor © as: 



Zn = 



N 



N 

E Km 

M=0 



y M(N — M) 



(5) 



Here (^) = 7V!/[M!(JV- M)!] is the binomial coefficient 
and z = exp (—A/TV) is the effective fugacity. One can see 
that Zjq is a polynomial in z, that is related to the control 
parameter A of the model, but z is not a size-independent 
quantity, and depends on the number of particles TV. 

Nevertheless we embark on studying zeroes of the EPF 
Zn- As a first step we find zeroes of Eq. J^J, considering z 
as a complex TV-independent variable. The results of our 
numerical calculations, using Mathematica, for three 
values of N are represented in Fig. Note that the 
order of the polynomial of Eq. J5j increases rapidly like 
TV 2 /4, and therefore we were not able to perform precise 
calculations of the roots beyond N = 71. 

With increasing TV these roots approach the unit cir- 
cle. One can argue that this should indeed be the case. 
First, let us associate with the EPF J5J the nonequilib- 
rium complex free energy density fN(z) — (1/TV) ln(Zjv). 
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FIG. 2: Zeroes of Eq. © in the complex A-plane, nearby 
the critical value A c = 2, for three values of N. The inset 
shows more roots for TV = 71. The continuous line is the 
analytical perturbative estimation 1121 of the line of zeroes in 
the thermodynamic limit, see main text. 



For large N and \z\ > 1 the EPF is dominated by the cen- 
tral term M = N/2 and hence fiq(\z\ > 1) ~ (N/4) lnz. 
Thus we define / W (z) = (iV/4) In z as a free energy in the 
\z\ > 1 region. On the other hand, for \z\ < 1 the dom- 
inant contribution are coming only from the M = and 
M = N terms, and thus /jv(M < 1) ~ 1/iV which in the 
thermodynamic limit defines f^ 2 \z) = for \z\ < 1. The 
above analysis indicates that the model for which Eq. (JjjJ 
is the EPF with z as a control parameter undergoes an 
abrupt transition at \z\ = 1. In the thermodynamic limit 
N — > oo, /at exhibits an infinite jump at the transition. 
To obtain the location of the zeros of Zjj in this limit we 
have to equate real parts of the nonequilibrium complex 
free energies on both sides of the transition ^(| . Namely, 
we require that 



Re/«(«) = Re/W(z) 



(2), 



(6) 



Using a polar representation z = re 1 ^ we immediately 
obtain that the only way to satisfy Eq. © is to have 
r = 1. Hence, asymptotically, the zeroes should be lo- 
cated on the unit circle, as confirmed by our numerical 
calculations. 

However, the model with control parameter 

(which has a transition with a jump of the effective free 
energy density) is quite different from the original urn 
model with A as a control parameter (which has a contin- 
uous phase transition). Therefore, in order to infer some 
information about the phase transition in the urn model 
we have to analyze the behavior of zeros of Eq. (J5J in the 
complex A-plane, that can be obtained from the zeroes 
in the z-plane using the relation A = — /Vln(z). Trans- 
formation of zeros into the complex A-plane is shown in 
Fig.0 

With increasing N the zeroes approach the critical 



point A c — 2 with a slope that is close to 7r/4, and with 
a vanishing density of zeroes. These numerical obser- 
vations seem to confirm the second-order nature of the 
phase transition 0. 

In the following we establish analytically the form of 
the line of zeroes close to the critical point, and the final 
result is presented as a continuous line in Fig. |21 As 
already mentioned, in the thermodynamic limit the EPF 
is dominated by the stationary state: 

(l/2-e)exp[-A(l/2-e)] 

= (l/2 + e)exp[-A(l/2 + e)]. (7) 

Below the critical point the leading term of for N — > 
oo is given by the central peak M — N/2: 
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On the other hand, for A > A c there are two leading 
contributions to Zn coming, respectively, from M — 
N(l/2-e) and M = N(l/2 + e), where e is the solution 
of the macroscopic stationarity condition J7J|. Therefore: 
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A > A c = 2. (9) 



Correspondingly, the effective free energy density / = 
limAr^ 00 (l/7V) ln(Zjv) H3 associated with this EPF is 

= —A/A + In 2, A < A c = 2, (10a) 
/ (2) = -.4(1/4 -e 2 )- (1/2 - £ )m(l/2-e) 

-(1/2 + e) ln(l/2 + e), A > A c = 2. (10b) 

Let us now consider the behavior of the EPF and of the 
effective free energy density as a function of the complex 
parameter A. Then the condition @, Re = Ref( 2 \ 
together with Eq. J7|) determine the line of zeroes in the 
complex A-plane. Note that now e is a complex variable 
obtained from the steady-state Eq. @. 

However, Eq. © is now too complicated to allow for a 
complete analysis of the zeroes line in the entire >l-plane. 
But we are mainly interested in the behavior of this line 
in the vicinity of the critical point A c — 2. Therefore, we 
shall look for a perturbative solution of Eqs. 10 and Q 
in the small real parameter a = Re A — 2 around A c . An 
inspection of these equations shows that Im A scales like 
a, while both Ree and Ime scale like a 1 / 2 . We are thus 
led to consider the following developments: 

A = (2 + a) + ia(ciQ + aa\ + a 2 a2 + ■ ■ .), (Ha) 
e = a 1 / 2 (x$ + ax\ + a 2 X2 + ■ ■ ■) 

+ia 1 / 2 (yo + ay 1 + a 2 y 2 + ■■■)■ (lib) 

We substitute these expressions in Eqs. © and Q, 
then solve them order by order in a. This leads to the 
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FIG. 3: The minimum distance between the zeroes of Zn and 
Ac — 2 in the complex A-plane as a function of N~ 1/2 . The 
continuous line is a least-square fit of the form a + &7V -1 / 2 (l + 
cTV" 1 ), where a, b, and c are fitting parameters. Note that 
the correction term N~ 3 ^ 2 is used in view of the small values 
of N that are accessible to calculations. Extrapolating to the 
N — > oo limit we obtain a/A c ~ 1%, i.e., very close to zero, 
that confirms the theoretical value A c = 2. 



following parametric expression for the line of zeroes in 
the j4-plane, in the vicinity of A c = 2: 



A = (2 + a) + ia [l + 0.6a + 0(a 2 )] 



(12) 



The result of this perturbative calculation up to 0(a 5 ) is 
represented by the continuous line in Fig. [21 Note that, 
indeed, the slope of this curve at the critical point is 
7r/4, sign of a second-order phase transition. Moreover, 
one can compute the density of zeroes on this curve in the 
vicinity of the critical point using the relationship |l0| : 
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where fj,(s) is the density of zeroes at a distance s from 
the transition point, distance measured along the line of 
zeroes. This gives in the vicinity of the transition point 
fi(s) = as + C(s 2 ), with s = V2a and a = 0.0045(1). 
The density of zeroes vanishes as a power law towards 
the transition point on the real axis - i.e., we recovered 
yet another characteristic of the equilibrium theory for 
second-order phase transitions. 

Up to now, we have shown that the zeroes of the EPF 
indeed provide information about the location and the 
type of the phase transition. But, as pointed out || cer- 
tain critical exponents are also encoded in the behavior 
of these zeroes. A similar conclusion can be drawn in 
our model. As shown in Fig. the distance \A — A c \ 
between the closest root A to A c = 2 in the complex 
A-plane decreases like Y -1 / 2 . A simple scaling argu- 
ment shows |5j that \A — A c \ should scale with the system 



size as N^ 1 ^ , where v is the correlation length critical 
exponent. However, our urn model is structureless and 
the correlation length does not seem to be a well-defined 
quantity. Nevertheless, we can implement a definition of 
the critical exponent v that is based on the finite-size 
scaling of moments of the order parameter |llj |. Since 
the probability distribution at the critical point for urn 
models is known || , using the standard prescription 
we obtain v = 2. Such a value agrees with the finite-size 
scaling observed in Fig. |3| 

IV. Conclusions. — We considered a simple stochastic 
model with state- dependent transition rates for which the 
EPF can be computed analytically and that exhibits a 
nonequilibrium second order phase transition. Our aim 
was to show that, although it is not a straightforward 
task to apply the concepts of the Yang-Lee theory, it is a 
remarkable and non-trivial fact that they still apply when 
the microscopic transition rates of the model are state- 
dependent. Indeed, all the other models that have been 
studied so far allowed for a rather simple transposition 
of the equilibrium theory. This was essentially due to 
the existence of some state-independent transition rates, 
and the possibility to express the EPF as a polynomial 
of these parameters. Our work therefore opens new per- 
spectives for the study of other nonequilibrium models 
with phase transitions, indicating that Yang-Lee's ideas 
apply to a much wider range of phenomena than those 
considered up to now. 
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